Thermodynamic phase transitions for Pomeau-Manneville maps 
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We study phase transitions in the thermodynamic description of Pomeau-Manneville intermittent 
maps from the point of view of infinite ergodic theory, which deals with diverging measure dynamical 
systems. For such systems, we use a distributional limit theorem to provide both a powerful tool 
for calculating thermodynamic potentials as also an understanding of the dynamic characteristics at 
each instability phase. In particular, topological pressure and Renyi entropy are calculated exactly 
for such systems. Finally, we show the connection of the distributional limit theorem with non- 
Gaussian fluctuations of the algorithmic complexity proposed by Gaspard and Wang [Proc. Natl. 
Acad. Sci. U.S.A. 85, 4591 (1988)]. 
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I. INTRODUCTION 

The major goal of statistical mechanics is to explain 
the macroscopic properties of complex systems in terms 
of a very small number of parameters by using probabilis- 
tic approaches. As is well known, its primary motivation 
was the study of thermodynamical properties of matter 
based on random behavior of their very large number of 
constituents (atom and molecules). Almost a century af- 
ter Boltzmann's seminal work [lj, such approaches were 
extended to dynamical systems theory (2|-|4J, a branch 
that is currently known as thermodynamic formalism Q ■ 
In this scenario, the chaotic dynamics of an ensemble of 
trajectories plays the role of randomness in the many- 
body dynamics, even for one degree-of-freedom dynami- 
cal systems. The thermodynamic approach has proven to 
be a powerful tool in the ergodic theory of hyperbolic and 
expanding dynamical systems [4J. Later, there has been 
growing interest, mostly by theoretical physicists, in ex- 
tending this approach to more general dynamical systems 
(see Q and references, therein), particularly those that 
exhibit fractal sets or some kind of intermittent 

behavior [HHI1]- 

Here we deal with phase transitions for Pomeau- 
Manneville (PM) maps x t +i — f{xt) where / takes the 
form 



f(x) = x(l + ax 1,a ) modi, 



(1) 



with a > and a > [17]. The remarkable character- 
istic of such systems is the intermittent behavior due to 
the presence of the indifferent fixed point x = 0, i.e., 
/(0) = and /'(0) = 1. It is important to stress that 
the global form of / far from x — is less relevant here. 
For example, systems behaving like ([]} on [0, where 
f(x*) = 1, exhibit the same statistical behavior of ([1]) 
since the map on [a;*, 1] is given by some well-behaved 
function /i such that f\{x„) — and /i(l) = 1. Systems 



of the type ([T]) have diverging invariant measure fi(x) 
near their indifferent fixed points for < a < 1. More 
specifically, the invariant density uj(x) of map ((T|), where 
d/j(x) = uj(x)dx, behaves as 



lj(x) ~ bx- 1/a , 



(2) 
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near x — 18]. Therefore, diverging measure regime of 
^ leads to a very slow laminar phase near x = al- 
ternating with fast turbulent one elsewhere. Due to this 
peculiarity, the dynamics of the system ([1]) exhibit subex- 
ponential instability of the type \Sx t \ ~ |<5x | exp(A a t") 
for < a < 1 [l9| . On the other hand, a > 1 leads to the 
finiteness of invariant measure, which is naturally related 
to the usual chaos and ordinary Lyapunov exponents. 

It is important to point out here the connection be- 
tween subexponential instability and the so-called "spo- 
radic randomness," a phenomenon initially studied by 
Gaspard and Wang [20] . These authors conjectured that 
the Kolmogorov-Chaitin algorithmic complexity C* for 
map (JTJ) is proportional to the number of entrances N t 
into a given phase space cell after a large number of it- 
erations. In this assumption, recently confirmed in (l9j 
by means of a Pesin-type indentity, the statistics of N t is 
ruled by non-Gaussian fluctuations involving Feller's re- 
newal results [2lT ] . Subsequently, thermodynamic phase 
transitions of PM map U]) for < a < 2 was studied by 
Wang [ll| employing the same approach of [2(| ■ It is also 
interesting to note that the sporadic randomness has not 
only been verified in PM intermittent maps (e.g., |19|), 
but also suggested as a distinguishing feature in weather 
systems [22|], noncoding DNA sequences [23j |. and some 
linguistic texts (24|. 

The purpose of this work is twofold. First, we will 
revisit the pioneering results of (l2{ . but now from the 
point of view of infinite ergodic theory [25] (for some 
applications, see also [lfj, [26|). In this first part some re- 
sults involving phase transitions of the so-called topolog- 
ical pressure are considerably improved. The topological 
pressure can be interpreted as a free energy density asso- 
ciated with the ensemble of trajectories. We also discuss 
the phase transition related to the Renyi entropy, extend- 
ing the results observed in [ill ] to the diverging measure 
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(nonergodic) regime of PM map (Q}. Finally, we show 
the connection between Feller's sporadic statistics and 
the infinite ergodic theory. Second, it aims at under- 
standing the phase transition problem from a dynamical 
point of view since singular behavior of thermodynamical 
quantities does not tell everything about dynamic char- 
acteristics of a system. The approach employed here also 
show us precisely what happens at each phase, particu- 
larly in the subexponential regime of map ([I]). 



II. TOPOLOGICAL PRESSURE 

In the thermodynamic formalism, systems of the type 
([1]) exhibit continuous phase transition, a situation where 
thermodynamic quantities vary continuously but not an- 
alytically when some external parameter of the system is 
changed. The paradigmatic example in the usual statisti- 
cal mechanics is the ferromagnetic material at zero exter- 
nal magnetic field: Ferromagnets lose their spontaneous 
magnetization when heated above a specific critical tem- 
perature T c and the derivative of the magnetization with 
respect to magnetic field (susceptibility) diverges at T c 
and zero field. As shown in 12], the many-body model 
that most closely resembles (fljis the Fisher-Felderhof 
droplet model of condensation [27j • 

Let us first consider the topological pressure P(/3), a 
kind of negative Helmholtz free energy of thermodynamic 
formalism, defined as @ 



PCS) = lim ylnZ t (/3), 



(3) 



with the corresponding partition function Z t (/3) given by 



Ztifi) = $>xp 



{**} 



t-i 



-^\n\f'(f k ( Xi ))\ 



fc=0 



(4) 



The set of points {xi} in Eq. is chosen as follows. 
First, consider a partition of phase space into disjoint 
boxes Ai so that transitions between nearest-neighbor 
configurations (£, i') are possible, i.e., /(Ai) n A# ^ 0. 
For each allowed sequence io, . . . , it-i, there is a subset 
A x (io, • ■ • , it-l) of phase space defined by 



A x (i , 



) = {x:f k (x)tA lk ,k = 0,...,t-l}. 



(5) 



The size of subsets Aj fc goes to zero as t — > oo. Then, for 
very large but still finite t, we pick a representative point 
Xi, one from each subset, and collect them as the set {x^}. 
It is important to stress, however, that the analytical de- 
termination of this set is not usually a practical task. We 
can circumvent this problem by replacing the summation 
over Xi by the integration over the conditional measure 
as follows (h is an arbitrary function) 



h(xi 



xe[o,i] 



da(x)\A x (i ,...,i t _i)\h(x), (6) 



where u(x) represents the measure of initial condition x 
and |A| denotes the Lebesgue measure of A. As t — > oo 
we have the following property [HI 

\A x (i , • • • i = l A *(*i' ■ • • »<t-i)ll/'(aOI. ( 7 ) 

where, on the right side of (JTJ), we can replace i\ by «2 
and by |/'[/(x)]||/'(x)|, and so forth, leading to 



t-i 



\A x (i ,...,i t - 1 )\ = l[\f[f k (x)} 



(8) 



fc=0 



Finally, we can rewrite Eq. @ by means of Eqs. and 
dSJ yielding 



da(x) exp 



t-i 

•E 

fc=0 



(1-/3)V In |/'[/ fe (x) 



(9) 



Before attempting to estimate Eq. © wc will make use 
of the Aaronson-Darling-Kac ( ADK) theorem [25[ , which 
is precisely applicable to PM systems of type ((J). For 
such systems, this theorem ensures that, for a positive 
function i? integrable over \i and an arbitrary measure a 
of initial conditions absolutely continuous with respect 
to the Lebesgue measure, we have 



fe=0 



vy c ' 



.(*)/ 



(10) 



as t — > oo, where ^ 7 is a non- negative Mittag-Lcffler 
random variable of index 7 € (0, 1] and expected value 
-E(£ 7 ) = 1. The corresponding Mittag-Leffler probability 
density function p 7 (^) is given by (2ll [20l] 



Ar(0 



r^xi+j) 



-9 J 



r 1 ^(i + 7 ) 

^1/7 



(ii) 



where g 1 stands for the one-sided Levy stable density, 
whose Laplace transform is g(u) = exp(— u<) (see (29l.l3G| 
for a detailed discussion) . For PM maps of the type (fTJ) , 
the index 7 is 



a, < a < 1, 
1, a>l, 



(12) 



whereas the coefficient c 7 (i) in Eq. pU|) takes the asymp- 
totic form 31 1 



C7O) 



1 /a\° sin(7ra) 

ba \a/ , ir a 
(blnty 1 , 

1, 



< a < 1, 

a = 1, 
a > 1, 



(13) 



as t — ¥ 00, recalling that b = lim a ^Q x 1 / a ijj{x). For a > 1 
we have introduced the Birkhoff ergodic case 7 = 1, for 
which the corresponding Mittag-Lcffler density reduces 
to pi(£) = 8(1 — £), as in the a = 1 case. Evidently, we 
can choose d = In \ f'\ in the ADK formula (|10l) . 
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Consider now the algorithmic complexity Ct of PM 
map (JU, valid for all a > 19]: 



C t (a0~5>|/'[/*(z)]|, 



(14) 



fe=0 



as < — > oo . Equations (fTU|) and (| 13[) lead to (see also [19j ) 
C t d 



(Ct 



^7 j 



(15) 



as t —> oo, where the ADK average (Ct) = /i 7 i 7 is given 
in terms of the average of generalized Kolmogorov-Sinai 
entropy [l!| 



fa/ry C-^ 



d»\n\f\ 



(16) 



Going back to Eq. (|9]), we can overcome the integra- 
tion problem over arbitrary a considering it absolutely 
continuous with respect to the Lebesgue measure. Such 
condition is sufficiently broad to assure that our results 
involving phase transitions typically do not depend on 
the initial condition distributions. This is somewhat sur- 
prising in the case of nonergodic regimes, i.e., < a < 1 
in the present case. After applying Eqs. (ITBl and (fT6|) in 
the ADK formula (TIT)]) , we have 



^P 7 (£)exp[-(/3-l)M 7 e]- (17) 



Note that Eq. (jTTJ) is just the Laplace transform of 
p 1 , which is given by the Mittag-Leffler special function 
E 1 (u), namely [H 



M u ) = E i( u ) = Yl 



[r(i + 7 H T 

n-y) 



^ r(i 

n=0 v 



(18) 



with u = (1 — /3)/i 7 i 7 . Now, considering the asymptotes 
E 7 (u) ~ 7 _1 exp(w 1 / 7 ) as u — > oo [HJ and E 7 (u) ~ as 
u — > — oo [33j , we have finally for all a > and /3 near 1 

P(/3) .|Mi-^^<i ; (19) 

observing that h 7 = (c 7 — > 0) for a = 1. Note that 
Eq. (|19l) is in accordance with the results first obtained 
in J12| for < a < 2 and later extended for all a > in 
[1 61 ] . It is noteworthy here that, unlike these approaches, 
the prefactor /i 7 in Eq. (|19p is obtained exactly, given by 
Eq. OH) for all a > 0. 



III. RENYI ENTROPY 

Let us consider now the phase transition related to the 
Renyi entropy [5j: 



K(0) 



1 



„ lim 

1-/3 t-Hx> 



(20) 



where = Pj(io, ■ ■ ■ )H-i) usually denotes the prob- 
ability that a randomly chosen initial condition (cr dis- 
tributed) on the phase space falls into Aj at time t— 1. In 
view of the fact that we are also dealing with nonergodic 
regimes (0 < a < 1), we set p^p such that 



E 



-1 a- 

k=0 T k 



(21) 



where pf\q = 1) = pf\ In Eq. ([21]) we consider the 
amount of time Tj spent in state Aj instead of its length 
\Aj\, which is usually considered for ergodic systems (for 
which Tj oc |Aj |). Then the partition function Z t takes 
the asymptotic form 



Z t (q) ~ exp[iP(g)] 



t-i 

k=0 



(22) 



Recalling that Y,j pf\o) ~ ex P[(! - l3)K(/3,q)t], we 
have for 5=1 



(23) 



also valid for ergodic systems @. From the topological 
pressure (fT9"f we then have 



^/ 7 (l-/3)- 1+1 / 7 , /3< 1, 



/?> 1. 



(24) 



Note that for 7 = 1, i.e., ergodic regimes, K(f3) — hies 
for j3 < 1, where ft-i = /i^g is the Kolmogorov-Sinai 
entropy, whereas K(j3) = for (3 > 1. Therefore, Eq. 
(l24l) extends for nonergodic regimes (0 < a < 1) the 
nonanalytic behavior of Renyi entropy at /3 = 1 observed 
in El. 



IV. ALGORITHMIC COMPLEXITY SATISFIES 
THE ADK THEOREM 

In [TJj], as well as in [2FJ], the algorithmic complexity 
Ct of a piecewise version of the PM map was considered 
as the random number of entrances N t into a given phase 
space cell (Aq) during t iterations of the map, i.e., Ct ~ 
Nt ■ The statistics p a of Nt employed is well known from 
Feller's renewal theorems [2L] , and it was applied in the 
estimation of P{/3). The accordance with Eq. (IT91 can 
be understood by observing that p a is, in fact, a Mittag- 
Leffler probability density function. The statistics of Nt 
for the case < a < 1 is given by [2l| 



(25) 



as t — > 00, where P Q and G a stand for the cumulative 
distribution functions of p a and g a , respectively. Apply- 
ing the change of variable q — r^ x / a , with r a = aT(a) 
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[29j, and after introducing the normalized random vari- 
able f = N t / (N t ), we have 

p a (N t )dN t ~ P*(t)dt, (26) 

as t — > oo, where (Nt) = c\t a /aT(a) and p a is the 
Mittag-Leffler density (fTTI) . For the case a > 1, but 
different from 2, we have [21| 

^ (jV t > c 2 t - cst 1 /^ ~ G K (q), (27) 

where k = a for 1 < a < 2 and n = 2 for a > 2. Now 
G K stands for the cumulative distribution function of the 
two-sided stable density g K [l2|. We can consider the 
same normalized variable £ of Eq. (|26l) . but now with 
(N t ) =c 2 t. Then Eq. ([27]) becomes 

Pa (Nt)dN t ~ i ff „ (^7^) (28) 

where e = (03/ 'c2)i _1+1/ goes to as i — > 00. This 
leads to 6(1 — £)e?£ on the right hand side of Eq. (l28l) . 
Already in the Gaussian case k = 2, e 2 is proportional to 
the variance, also leading Eq. (|28|) to the same Dirac 6 
function. The same occurs for a — 2, where we also have 
Eq. (|28p for k = 2, but replacing C3 by c^Vhii and also 
leading to e — > as i — > 00. 

V. FINAL REMARKS 

We revisit here the problem of thermodynamic phase 
transitions for PM maps ((TJ) by using the infinite crgodic 
theory, in particular the ADK theorem. The topological 
pressure P(/3) and Renyi entropy K(/3) are calculated ex- 
actly for such systems, exhibiting both phase transitions 



at the same critical value j3 c = 1. Our results also shed 
some light on the role of the measure of initial conditions 
(7 in the calculation of these thermodynamic functions. 
Such quantities are invariant under a since it is abso- 
lutely continuous with respect to the Lebesgue measure. 
This result is somewhat surprising in the case of the non- 
ergodic regime of the PM map ([T]), showing once more 
the strength of the ADK theorem. 

From a dynamical point of view, the thermodynamic 
formalism allows us to obtain important quantities that 
characterize nonlinear systems. We can mention, for in- 
stance, the Pesin formula relating Lyapunov exponent 
A to the Kolmogorov-Sinai entropy hi = hxsi namely 
h K S = [P(P) ~ P'W)]^i- = A| For nonergodic 
regimes < a < 1, however, the topological pressure (fT9)l 
gives us the trivial relation Iiks — A = 0. In fact, the dy- 
namic instability is stretched exponential for such cases, 
of the form \Sx t \ ~ \Sxq\ exp(X a t a ) , rather than exponen- 
tial a = 1. For such cases we can consider, once more, 
the infinite ergodic theory approach. It has recently been 
shown that the Pesin relation can be extended in a non- 
trivial way provided one introduces a convenient subex- 
ponential generalization of the Lyapunov exponent and 
Kolmogorov-Sinai entropy Moreover, the general- 

izations of such quantities behave like Mittag-Leffler ran- 
dom variables with h a as the first moment [lj| [26[ . A 
quest for new constitutive relations involving P(f3) that 
lead directly to these results in the thermodynamic for- 
malism probably deserve further investigations. 
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